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Abstract 

The mechanism by which a threshold may capture a resonance is examined. It involves a 
threshold cusp interfering constructively with either or both (i) a resonance produced via 
confinement, (ii) attractive t- and -u-channel exchanges. The /o(980), X(3872) and Z(4430) 
are studied in detail. The /o(980) provides a valuable model of the locking mechanism. 
The X(3872) is too narrow to be fitted by a cusp, and requires either a resonance or virtual 
state. The Z(4430) can be fitted as a resonance but also can be fitted successfully by a 
' cusp with no nearby resonant pole. 

^ i PACS numbers: 12.39.Mk, 13.25.Gv, 13.25.Hw, 14.40.Lb. 14.40.Nd 

1 Introduction 



It is well known that several meson resonances appear at or close to thresholds. Examples are 
/o(980) and a (980) at the KK threshold, / 2 (1565) at the uju threshold and A (1430) close to 
Ph| the Kr]' threshold. In all these cases, decays to mesons are via S-waves. Further possible ex- 
amples are A(3872) at the D(1865)£>*(2007) threshold, Z(4430) close to the £>* (2007)1)1 (2420) 
threshold and F(4260) close to the -D(1865).Di(2420) threshold. [The DD charge conjugate corn- 
er . bination is tacitly included throughout this paper]. Amongst baryons, Pn(1710) and Pi3(1720) 
appear close to the Nuo threshold (1720 MeV), and A c (2940) appears close to the L>*(2007)A 
(*C) \ threshold. 

All these cases have been discussed as 'molecules', bound a few MeV below the thresholds. 
There is much argument whether they should be viewed as bound states of mesons or as quark 
^ ■ configurations, e.g. ccnn in the case of A (3872). 

The first objective of this paper is to draw attention to a known but frequently overlooked 
mechanism which attracts resonances to thresholds. The /o(980) is studied in depth as one 
specimen. The second objective is to point out the role of zero-point energy. This provides 
^ ■ some distinction between mesonic molecules (uncoloured) and quark configurations, which are 
coloured and therefore confined. Thirdly, an improved but convenient form of the Flatte formula 
PP is proposed for the case of sharp thresholds like KK and other cases cited above. 

Consider as an example /o(980) and its decay to KK. The conventional Flatte denominator 
for this resonances is 

D(s) = M 2 -s-tJ2^ (1) 

i 

where M is resonance mass, Imlli = gfp;, gf are coupling constants to decay channels and pi 
are phase space factors. However, further dispersive terms Rellj(s) are required in D(s): 

-ReIIM = -lp rd^WflW, (2) 
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where are it and K masses and P denotes the principal value integral. In fact, the terms 
igfpi of Eq. (1) arise from the pole at s' — s in Eq. (2). 

Fig. 1 shows Hkk(s) and Q 2 kkPkk{s) near the KK threshold, using a form factor e~ 3k , where 
k is centre of mass KK momentum in GeV/c. There is a cusp in U(s)kk at the threshold. The 
fact that it is positive definite at theshold signifies additional attraction appearing there. A 
minor technicality is that IT(s) goes negative after the peak in p(s), since the dispersion integral 
is close to the gradient of p(s). Eventually it returns slowly to zero for large s. 
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Figure 1: Ren^x(s) and 9kkPkk(s) for / (980), normalised to 1 at the peak of g 2 KK PKK- 

If the cusp is superimposed on attraction from another source, for example meson exchanges, 
a resonance can be generated by the peak in Re II; if the attraction is not quite sufficient to 
produce a resonance, there is a virtual state instead. The <r, k, a (980) and / (980) make good 
candidates for such molecules. Otherwise, all known cases except Z(4430) can conservatively 
be accomodated as regular quark resonances whose masses have been perturbed to synchronise 
with thresholds. 

Near threshold, the wave function of the resonance has a long tail, like the deuteron. The 
long-range tail is purely mesonic; any contribution to the wave function from coloured quarks 
is confined in the resonance at short range or in the decay mesons. In the confined part of the 
wave function, quarks have kinetic energy k 2 /2m = — (/i 2 /2m)V 2 \l/; here k and m are quark 
momentum and effective mass and ^ is the wave function. If a bound state is to be created, 
extra potential energy is needed to compensate the kinetic energy of the confined particles. This 
again pulls the resonance towards the threshold. 

Jaffe proposed that the / (980) and a (980) are members of a 4-quark nonet with composition 
ss(uu ± dd) [2]. This may be correct, but a long-range tail of KK is unavoidable. 

Section 2 reviews the formalism. This leans heavily on a perceptive paper of Tornqvist [3], 
written in 1995, before many of the examples cited above were known or accurately measured. 
For the case of Z(4430), some changes need to be introduced into his formalism. 

Section 3 examines the /o(980) in detail. Perturbations around its fitted parameters illustrate 
the general way in which a threshold captures a resonance. Section 4 examines X(3872), which 
peaks within 0.6 MeV of the DD* threshold. It could be a mesonic molecule, as Tornqvist 
[I], Close and Page [5] and Swanson [6] have argued. Its very close relation to the threshold 
points strongly to the cusp mechanism playing a decisive role. A likely possibility is that the 
cc 3 Pi radial excitation has been captured by the DD* threshold. If that state were to appear 
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elsewhere, the alternative explanation is a mesonic molecule or virtual state. Present data simply 
require a second-sheet pole. 

The Z + (4430) reported by Belle [7j appears as a peak in the mass spectrum in B — > 

ty'ir K. The remarkable feature of this peak is that it has isospin 1, unlike regular cc states, 
so it is exotic. Belle fit it as a resonance with mass M = 4433 ± 4 ± 1 MeV and width 
T = 44ljgtii MeV. Rosner [S] and also Meng and Chao [§] point out that its mass is close to 
that of D*(2010) + _Di(2420). Maiani et al. [10] interpret the peak as a diquark-antidiquark 
state. For an S-wave D*D\ combination, its spin-parity J p is 0~, 1~ or 2 _ . It is discussed in 
detail in Section 5. It is easily fitted as a resonance. However, it can also be fitted as a cusp 
due to strong de-excitation of D*(2010)/)i(2420) to lower configurations of DD and D*D. 

Section 6 makes remarks about a number of the other states associated with thresholds, so 
as to fill in the broader picture. Section 7 summarises conclusions. 



2 How to parametrise a cusp or threshold resonance 

Tornqvist's formalism [5] will be followed here with some modifications. He considers the possi- 
bility of overlapping resonances such as /o(980) and /o(1370); that is a detail not required here, 
simplifying the notation. His formula for the amplitude connecting channels i and j (e.g. tttt 
and KK) is then 

A tJ (s) = T u (s)^ = Gl(s)PG 3 (s), (3) 

where P is the propagator; T is the usual T-matrix with phase space factored out. The Gi are 
related to coupling constants via 

G 2 (s)=g 2 Pl (s)F 2 (s)6(s-s th ,) } (4) 

with s t h,i the threshold for channel i; Fi(s) is a form factor, taken here as exp(— k 2 R 2 /6) with 
R = 0.6 fm. Then F 2 ~ exp(— 3A; 2 ) with k in GeV/c. For narrow resonances, the precise 
s-dependence of F 2 is not decisive but begins to matter for broader structures like Z(4430). 
For a resonance, 

P-\s) = M 2 - s-U(s). (5) 

The definition of n has the opposite sign to that used by Tornqvist, to avoid a multiplicity of 
minus signs in formulae which follow. The unitarity relation Im A = AA^ gives 

ImH(s) = 5:g 1 2 p i (s)Ff(s)^s-s th , i ). (6) 

i 

Using analyticity, 

Ren(s) = ip f" ds-I^. (7) 

7T Js thA S — S 

The term inside the summation of Eq. (6) is positive definite , leading always to attraction at 
and below threshold, though the dispersion integral can change sign above threshold. 
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2.1 Useful resonance formulae 



If F 2 (s) = 1, the integral in (7) diverges. However, with two subtractions, it can be evaluated 
analytically with the result [12] 



ReIT(s)i . pi I- pi 

2 — = Ji = — ln T7 ' s ^ S th,s 

9i 7T 1 + pi 



where 



'4m? — s 2vi , , . 
tan Vi, s<s thi (9) 

S 71 



Pi = 2ki/y/s, S>S th ,i (10) 
Vi = 2\ki\/y/s, s < s th ,i (11) 

and hi are momenta of decay particles in the rest frame of the resonance. The first term of Eq. 
(9) is the usual Flatte extrapolation below threshold. The contributions to ji near threshold 
from remaining terms are 

nji = -2p 2 l -(2/3)pt + ..., s>s thji (12) 
= -2v 2 + (2/3H 4 + . . . , s<s th>i . (13) 

The first term in each is symmetrical about threshold, but the next term leads to an asymmetry 
away from the threshold. It is also important to realise that the two subtractions used in 
arriving at this result allow IT to contain an additional linear dependence on s. A constant 
term is certainly required to make the integral of (7) positive at threshold. 

In illustrations given later, the principal value integral will be evaluated numerically. However, 
the algebraic forms (8) and (9) improve on the usual Flatte form for a resonance near threshold, 
but are equally convenient to use. A constant term can be absorbed into M 2 of the Breit-Wigner 
amplitude giving 

D(s)=M 2 -s-J2gK^ + ^)- (14) 

i 

For small v\ below the opening threshold, the second term of Eq. (9) contributes 

71 71 \ S J 

Note the explicit appearance of the factor (Am 2 K — s) in equation (15) and also in the first term 
of Eq. (9). The cusp mechanism is contributing to the resonance through a term like (M 2 — s) 
but peaking at the threshold. Also note that, when fitting data, there is a strong correlation 
between the term M 2 — s and —g\ji- If there are no data on the KK channel to determine g\ 
directly, data below the KK threshold give almost no determination of g\ in view of the poorly 
known form factor F|. 

The fit to BES data on J/f -> 07r + 7r~ and <pK + K~ [llj changes very little when the Flatte 
formula is replaced by Eq. (9). The (pKK data deterine g 2 well, and a modest change in M 2 
compensates for the change in the formula. 
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If there is a linear subtraction term in addition to the constant term, 



D(s) - M 2 -(l+/3) S -$> 2 (j 4 + w ) 



(16) 



M 2 



1 



Y.9 2 i{3i + ipi) 



(17) 



1 + 



— s — 



This leads to a renormalisation of the mass and width, as observed for /o(1370) in a recent re- 
analysis including dispersive terms arising from the opening of the 47r channel [13] . For practical 
purposes, the factor 1/(1 + /?) renormalises mass and width terms in (17). 

2.2 No resonance and approach to resonance 

If there is no resonance, it is convenient to replace M 2 — s of (14) by a constant M 2 , hence 
keeping dimensions unchanged. In Section 5, it will be shown that the K-matrix approach leads 
to further terms neglected by Tornqvist, but they can safely be ignored in the cases of f (980) 
and A (3872). Standard effective range theory replaces M 2 by M 2 — jk 2 . As 7 increases, a 
virtual state approaches the threshold and smoothly becomes a bound state - or resonance if 
there are open channels. 

The cusp mechanism gives maximum attraction at the threshold. A resonance close to thresh- 
old necessarily has a long-range tail, mesonic in character. Consider as an example mixing be- 
tween cc and a (cn)(cn) configuration, appropriate to the discussion of X(3872). The mixing is 
given by the familiar eigenvalue equation 



where H\\ and if 22 describe isolated cc and (cn) (cn) configurations and V describes mixing. This 
mixing pushes the cc state down in mass. This is the well known Variational Principle which 
minimises the eigenvalue when states mix. It is closely analogous to formation of a covalent 
bond in chemistry. The cn configuration is different in detail from an atom made of a proton 
and electron, but the principle is the same. The detailed dynamics of the decay of a confined 
state through the confining barrier is a key missing detail at present; so is the possible role of 
di-quark interactions in the 'molecule'. 

The meson-meson component has t or u channel poles, but slow variation with s. The linear 
combination with cc has minimum zero-point energy when the resonance is centred on the 
threshold and is narrow. Both the cusp attraction at threshold and the effect of zero-point 
energy provide feedback tending to lock the resonance to the threshold. In the next section, the 
/o(980) will be used as a model to examine this locking mechanism numerically. 

The dispersive effect at the threshold may be an unfamiliar effect. Vacuum polarisation arises 
in this way. An example from classical physics is a tsunami - a travelling cusp. As a tsunami 
approaches a beach, attraction into the wave-front drains water from the beach well in advance 
of the wave and gives advance warning of the approaching wave. 
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3 The / (980) 



Can we be sure that the /o(980) is really a resonance and not a virtual state or cusp? Cern- 
Munich data on tttt elastic scattering [14] were first fitted to a simple Breit-Wigner resonance of 
constant width and yielded a good fit with a phase increase of ~ 180° over the KK threshold. 
Crystal Barrel data on pp — > 3tt° are a further delicate source of information [15]. The /o(980) 
appears as a visible dip in the irn mass projection. These data provide accurate phase informa- 
tion on / (980) via interferences between the three 7r°7r° combinations all over the Dalitz plot. 
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Figure 2: (a) The Argand diagram for the tttt S-wave from [13], (b) enlarged comparison with 
the /o(980) mass region; crosses indicate errors for free fits to real and imaginary parts of the 
amplitude in 10 equal bins of s from 0.84 to 1.08 GeV 2 ; dotted lines indicate the movement of 
free fits from that with a Flatte formula; numbers indicate masses in GeV. 

Fig. 2(a) shows the Argand plot of the tttt S-wave amplitude from a recent re-analysis of the 
data [T3] using a Flatte formula for / (980). As a test of the resonance hypothesis, these data 
have been refitted allowing complete freedom for the real and imaginary parts of the /o(980) 
amplitude in 10 equal bins of s wn from 0.84 to 1.08 GeV 2 . Fig. 2(b) compares results (crosses 
indicating the errors) with the fit of [13] shown by the smooth curve. The errors allow for a 
mass resolution of ±4 MeV, which has quite a large effect near the KK threshold. 

It is clear from Fig. 2(b) that the phase shift of /o(980) increases by ~ 90° from ~ 900 MeV 
to the KK threshold at 991 MeV (the average mass for K + K~ and K°K°). This supports the 
resonance interpretation. 

As a further test, the fit has been repeated using Eqs. (8) and (9), which improve on the 
Flatte formula. Also, separate thresholds have been included for K + K~ and K°K°. There 
is a small improvement of 23 in log likelihood (statistically 5.5 standard deviations). There is 
only a small shift in the second sheet pole, which moves from 998.4 ± 4 — z(17.2 ± 4) MeV to 
(1003.9 — il6.5) MeV with identical errors. [For the second sheet pole, the signs of p2 and ji are 
reversed; for the third sheet pole, they are as in Eq. (14)]. Because the parametrisation of the 
cusp is now included, the distant third sheet pole moves from (851 — z418) MeV to (1171 — i592) 
MeV; this is because the cusp formula modifies the shape of the resonance away from the KK 
threshold. 

It is of interest to use the /o(980) as a model to examine the movement of the second-sheet 
pole as the parameters of the resonance are varied. This reveals the general features of the way 
the resonance forms. These general features are likely to be similar for other cases. 
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M g\ Pole 
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Table 1: Movement of the second sheet pole of /o(980) with g| using a pure cusp formula. 



A fundamental point is whether the cusp alone can produce a pole at or very close to the 
KK threshold. The answer is no. For this test, the term — s of the Breit-Wigner denominator, 
Eq. (14), is removed, again including KK mass differences. A second sheet pole does appear, 
but well below the KK threshold. Table 1 shows the pole position for a variety of values of g\. 
For small , the pole moves rapidly away from the KK threshold. The reason the pole always 
lies below the KK threshold is that it is closely related to the term 1\k^\j yfs of Eq. (11) which 
increases below the KK threshold with |fc 2 |. 

The conclusion is that formation of the resonance requires an additional source of attraction, 
e.g. meson exchanges. Janssen et al. [IB] were able to generate a resonance from K* and p 
exchanges between tctc and KK. Likewise, the a and k poles appear from meson exchanges in 
the calculations of Caprini et al. [T7J and Biittiker et al. [H] using the Roy equations. 

The attraction due to these long-range forces makes the real part of D(s) pass through 
at some mass. It is of interest to see how effective the cusp is in attracting the second sheet 
pole when this crossing point varies. Table 2 shows the pole position as M of the Breit-Wigner 
formula is varied. From the first few entries, one sees that the pole moves a long way if M is 
far from the KK threshold. The cusp is effective in attracting the resonance over a mass range 
much larger than the 7T7T width. From M = 0.9 to 1.10 GeV, the pole stays quite close to the 
KK threshold. For M outside the range 0.5 to 1.1 GeV, it disappears. The pole is more easily 
pulled up to the KK threshold than down to it. Over the range of values in the Table, the 
distant third sheet pole moves only by small amounts. 

The bottom entries in the Table show how the second sheet pole varies with g\. As it 
decreases, the width of the pole decreases more rapidly; as g\ — > 0, it becomes a bound state 
0.15 MeV below the K + K~ threshold (the lower of the two KK thresholds). This bound state 
survives unchanged for zero g\ for any value of M over the range 0.5 to 1.1 GeV and beyond. 

The general pattern which emerges is that a cusp superimposed on attraction due to either 
meson exchanges or quark confinement can lock the resonance close to the threshold for a wide 
range of parameters M and g\ . The cusp acts as a trigger for the resonance. This is particularly 
the case when g\ is small. The locking mechanism is obviously inhibited if there is repulsion in 
meson exchanges. 

Tornqvist [3] gives a formula for the KK component in the wave function of any resonance 
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M g\ Pole 
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0.02 
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Table 2: Movement of the second sheet pole of / (980) from Eqs. (8) and (9) as its parameters 
are varied. In all cases, g\ = 0.875 GeV 2 . 



close to a threshold. For the present case, it may be written: 

, \qqqq> +Eii(d/ds)ReIl i (s)] 1/2 \KK > 

l + E i (d/ds)ReU i (s) ' 1 1 

this was evaluated in Ref. [19] and gives ~ 60% KK component for / (980). That reference 
also shows that decay branching ratios favour an abnormally large KK component in /o(980). 

Model calculations of a similar nature have been made by van Beveren and Rupp [20]. They 
adopt a transition potential coupling confined states in a harmonic oscillator potential to out- 
going waves, with a matching at the transition radius ~ 0.65 fm. The a, k, a (980) and / (980) 
emerge from the continuum as their coupling constant to confined states is increased. The great 
merit of their model is that it is straightforward to follow the movement of poles as this coupling 
constant varies. They give graphic illustrations of the movement of the k and a o (980) poles with 
coupling constants. The movement of the a, k, a (980) and / (980) poles as a function of cou- 
pling constant is also tabulated in Ref. [21]. The model reproduces fairly well the amplitudes for 
all these states with a universal coupling constant. The do does not appear at the rjTc threshold 
because of the nearby Adler zero. The a, k and a® all become bound states if the coupling 
constant is increased by a factor 2.5. 

Jaffe proposes that mesons can be divided into 'ordinary mesons' which decouple from scat- 
tering channels as the number of colours iV c — > oo and 'extrordinary mesons' which disappear 
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in this limit [22] ■ He draws attention to the work of Pelaez on the N c dependence of the chiral 
Lagrangian and its predictions for low mass tttt scattering [23] . 

A final comment is that the phase shift for a goes through 90° very close to the KK threshold. 
This may not be accidental. The a has significant coupling to KK and it is possible that the 
cusp mechanism is sufficient to tie the 90° phase to the threshold as well as that of / (980). 



Fig. 3 shows j2 and p 2 for the DD* threshold. The result is similar to Fig. 1. Attempts to fit the 
data with a bare cusp (i.e. with M 2 — s of the denominator replaced with just M 2 ) can reproduce 
the very narrow width observed in decays to pJ/^ only with very fine tuning of parameters. 
There is then a second sheet pole almost at the DD* threshold. However this narrow width 
is incompatible with data on decays to DD* shown below in Fig. 4. The conclusion is that a 
resonance or virtual state is required 

The cusp can however capture a nearby cc 3 Pi state. The alternative possibility is a molecule 
generated by n exchange. 



The data with the best mass resolution are from Belle and Babar on decays to pJ/^ [21] [2S] 
and DD* [26] [21]. These have been fitted by Braaten and collaborators [29] and by Hanhart 
et al. [30]. They favour respectively a weakly bound state and a virtual state. Both omit the 
dispersive effect discussed here, but assume an attractive interaction which generates a nearby 
second sheet pole. 

It is necessary to fit the narrow width observed in pJ/^ simultaneously with the ~ 3.5 MeV 
mass difference between the peaks observed there and in DD*. The Belle data are shown in 
Fig. 4 with a variety of fits. For pJ/^ the improved data set shown by Olsen is used [31]. 

The mass distribution in pJ/^ (taken as channel 1) and DD* (channel 2) are: 
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Figure 3: j_dl>*( s )) Pdd*( s ) f° r A(3872); the vertical scale is arbitrary. 



ds 
dN 2 



ds 




(20) 



(21) 



D(s) 



M 2 -s-iM(T 1 +g 2 2 lp l (s)+j 2 /i]) 



(22) 
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Figure 4: Fits to Belle data for (a) and (c) X(3872) -> pJ/V, (b) and (d) DD*; the upper 
panels are fitted with a virtual state and the lower two with a resonance. 



for a resonance; for a virtual state, — s in D(s) is replaced by — 7&|. The parameter A is an 
overall normalisation constant. A constant Ti describes the width to channels other than DD*. 

The peak in DD* is naturally higher in mass than in pJ/^ because of the factor p2 in the 
numerator of Eq. (21). However the width of the J/ty peak cannot be too small if the difference 
in peak positions is to be accomodated. The branching ratio BR[pJ/^f]/BR[DD*] = 0.08±0.04 
including allowance for D* — > D^ as well as Dir. The branching ratio BR[uo J/ty]/ BR[pJ/ty] = 
1.0 ± 0.4 ± 0.3 [32]. A fit with these decay modes gives only a 1.5 MeV mass difference between 
pJ/'fy and DD*, but is sensitive to the assumed mass resolution; here a moving average is 
calculated over the 4.25 MeV bins quoted by Belle for DD* data. 

There may well be other unobserved decays, e.g. to [xict]l=o and/or [t]c<j}l=i, where L 
is the orbital angular momentum in the decay; they would be hard to detect experimentally. 
Contributions from these channels seem likely since decays to pJ/^ are isopin violating and 
the phase space for ujJ/^ is very limited. The splitting between pJ/^ and DD* peaks is 
easily increased to 3 MeV by allowing extra decay channels with a rate a factor 4-8 larger than 
pj/ty + uJ/ty. 

Figs. 5(a) and (b) show fits with a virtual state; Figs. 5(c) and (d) use a resonance. There is 
little to choose between them. Presently the only way of making a definite distinction would be 
to determine the phase directly from interferences with other components in Dalitz plots; the 
virtual state has a falling phase above the peak. A minor detail is that Fig. 4(d) show a second 
fit as a dashed curve, assuming there is some background in DD* at high masses. In view of 
the uncertainties, a full list of fitting parameters is not useful. For the full curve of Fig. 4(d), 
M = 3871.8 MeV, T = 8 MeV and equal branching ratios are used to DD* and the sum of 
other channels. The full-width at half- maximum is however only 4.3 MeV because of the cusp 
contribution to the resonance. 

Belle report a peak in DD* at 3942±£ ± 6 MeV with T = 37^ ± 8 MeV [28]. A spin- 
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parity analysis is important. It could be the 3 Pi radial excitation, in which case X(3872) is 
to be interpreted as a molecular state generated by the cusp mechanism and attractive meson 
exchange. Alternatives are J p = _ , 1~ or 2~, though these would be suppressed near the DD* 
threshold by an L = 1 centrifugal barrier for decay. 

5 Z(4430) 

The Z + (4430) observed recently by Belle is important as a clear candidate for an exotic meson. 
It appears at the D*(2010)-Di(2420) threshold within errors. It has a modest width ~ 44 MeV 
comparable with the width of D\ itself (25 ± 5) MeV. 

Unlike X(3872) it has many likely de-excitation processes 

. o- -> [ J D J D (2308)] L=0 , [DD*] L=1 and [D*D*] L=1 , 

• 1 -> [D*D (2308)} L=0 , [DD] L=U [DD*] L=1 and [D*D*] L=1 , 

• 2- -> [DD*] L=1 and [D*D*] L=1 . 

The essential point is that these de-excitation processes unavoidably lead to a strong threshold 
cusp unless meson exchanges are repulsive for all J p . 

It is necessary to fold the width of the _Di(2420) into the calculation of the cusp. Consider as 
an example decays of D*(2010).Di(2420) to D*(2010).D(1865) with L = 1. For fixed total s of the 
initial state and fixed mass s\ of Di, the intensity is proportional to (2fc 1 / v /i)S 1 (/c 1 i?)F 2 (/c 1 i?) 
(ignoring the very small widths of D and D*)\ here k\ is the centre of mass momentum in the 
final state and B\ = k\R 2 /{I + k\R 2 ) is the centrifugal barrier factor for L = 1 decay. Since k\ 
is large (~ 1070 MeV/c), B\ ~ 1. Including the phase space for the initial state and integrating 
over the line-shape of -Di(2420), the phase space factor for the whole process is 

p(s)= rfgi - — F (kiR), (23) 

J{M D +M D »y \D{sx)\ l v /ssT 

where k 2 is the momentum in the D*Di channel. This integral is easily done numerically. The 
important dependence is on k 2 , F 2 {k 2 R) and D(si). 

Fig. 5(a) show the result for p(s) and its dispersion integral g 2 j(s) with a subtraction at 3.95 
GeV, well below the D*D\ threshold. The full curves of Fig. 5 show p(s) (normalised to 1 at 
its peak) and its dispersion integral g 2 j(s) (normalised accordingly). The dashed curve shows 
the result of a large change in F 2 to exp(— 6k 2 ). 

The Belle data show a Z(4430) peak accurately coincident with the DD* threshold. The data 
can easily be reproduced as a resonance, and the fit is shown in Fig. 5(b). Close [22] argues in 
favour of binding by 7r exchange, though Liu et al. reach the opposite conclusion [34J. Li, Lii and 
Wang raise the interesting possibility of a nonet of strange relatives at the D*(2112)Z)i(2420) 
and D*D sl thresholds near 4540 MeV and the D*D sl threshold at 4647 MeV. 

The width of _Di(2420) reduces the sharpness of the cusp, though this could be more than 
compensated by the large number of decay channels into which D*D\ may de-excite. A cusp 
plus meson exchange could generate an exotic resonance, but the net attraction must overcome 
zero-point energy. We proceed to test whether an alternative fit is possible using a bare cusp. 
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Figure 5: (a) g 2 (D*D 1 ) and ReIl(D*D 1 ) calculated from de-excitation of D*D 1 to [D*D] L=1 . 
The dotted curve shows the effect of doubling the exponent of the form factor, (b) Fit to Belle 
data with a resonance; the dashed curve shows the predicted peak in D*Di, but scaled upwards 
for visibility 




5.1 A subtlety in the equations 

At this point, a digression is needed to discuss a missing element in Tornqvist's equations. For a 
2-channel system (which is sufficient to reveal the essential point), the usual expression relating 
the T-matrix and .fT-matrix elements is: 

T " = 5 ( t t ) ( (1 ":f 22 ' ff* ) < 24 > 

D \ K>2\ K 22 / \ — %K 2\ U — IKll ) 

B = 1 - iku - ik 2 2 + (&12&21 - &11&22), (25) 
where kij = ^Jkikjdij. Using time-reversal invariance, k\2 = &2i- Then 



Q22P2 ~ iPlP2(a 2 l2 - Q 11 Q 22) 

1 - ia n pi - ia 2 2P2 + £1^2(012 - 011022) 



u-22^2 - ^1^2V,"12 ~~ "11"'22; 

T 22P2 = - — — — : — — — : — , — f 2 — \ ( 26 ) 



T12VP1P2 = z : : 7—2 r, (27) 

1 - ianpi - ICL22P2 + P\P2\a>\2 - 011022) 

T11 is like T22 with indices interchanged. In order to arrive at a form T = iV/D where N is 
real, it is necessary to multiply top and bottom of T22 and Tn by the complex conjugate of the 
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numerator. Note that the resulting denominator D is then different from that of T 12 . So the 
result does not describe a resonance, for which D should be common to all channels. For all 
denominators to be the same, the requirement is that a\ 2 — 011022 = 0. This is satisfied at a 
pole. This term was omitted by Tornqvist and from the equations of Section 2 because there 
are indeed second sheet poles very close to the experimental peaks. In general however, this 
term is not zero. If it is completely dominant, Tupi and T22P2 —* i- 111 particular, this arises if 
T\2 is large. That is to be expected for de-excitation of D*D\ to many DD channels with lower 
thresholds. 

Fig. 6(b) shows the fit to Belle data with a bare cusp and Fig. 6(a) shows the Argand diagram 
for Ti 2 . Below threshold, the cusp effect creates attraction, pulling all amplitudes around the 
periphery of the unitary circle towards the threshold. At the threshold, the amplitude curls 
towards the centre of the Argand diagram. As the inelasticity grows, amplitudes for individual 
decay channels move back to the origin, though their total goes to i. For this fit, there is no 
second or third-sheet pole anywhere in the vicinity of the cusp. 
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Figure 6: Fit with a bare cusp: (a) Argand diagram, (b) fit to Belle data 



Some further detail is needed on exactly how data are fitted. The tails of a cusp depart 
somewhat from a resonance form. To allow for this, it is necessary to make small adjustments 
to the background fitted by Belle. This background is taken to be linear in s (slightly higher at 
low s than high s), multiplied by the phase space factors Q 1 for Z — > ip'n and Q2 for B KZ. 
The calculation of is done for Z — > [DD*]l =1 as a representative example. The dispersion 
integral is strictly over P1P2, but the dependence on p 1 is in practice small. It is assumed that 
ai2 a 2 2 though the coupling constant to DD* is fitted explicitly and so is an, which is taken 
as a constant since it is insensitive to form factors. The denominator of T12 for the Wir channel 
is the same as for all other channels. 

How can experiment resolve a cusp from a resonance? The difference in line-shape is small 
and easily confused with experimental backgrounds. The essential feature of a resonance is that 
the phase shift increases by 90° up to the centre of the cusp. The detailed s-dependence on 
Fig. 6(a) is different from a resonance, but the difference is delicate. This can in principle be 
measured from interference with other components in the Dalitz plot, but would require a large 
increase in statistics. 

The Z(4430) may be a resonance. The Argand diagram of Fig. 6(a) resembles a resonance. 
If meson exchanges are attractive, it is quite possible that they will combine with the cusp 
mechanism to generate a resonance. Mixing with a genuine diquark-antidiquark system is less 
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likely, since the long-range tail of the wave function must be purely molecular. 



5.2 Other points 

Rosner and also Meng and Chao remark on the experimental fact that no signal is seen for 
Z(4430) — > J/tyn despite larger phase space. This may be due to a cancellation in the matrix 
element for D*D\ — > J/^/n. The momenta for production of J/\I/ and \&' are large: 1130 and 
670 MeV/c respectively. The J/\l/ and ^' wave functions are multiplied by a factor exp(ifcr) 
which has nodes at 0.27 or 0.43 fm in the two cases. 

In more detail, the matrix element for Z — > ^'ir in the rest frame of the is 



using for the initial state the long-range part of its wave function, expressed in terms of spherical 
Bessel and Neumann functions. This is only a first approximation and the short range part of 
the wave function is unknown. The appearance of ri\ sin 5 signifies the fact that the radial wave 
function is sucked in by the attractive interaction due to the cusp, causing the phase shift S. 
For a strong cusp due to many open channels, n\ sin S should dominate. For the kinematics of 
J/ty production, n\ has a node at 0.38 fm, leading to a distinct cancellation within the matrix 
element. For production, n\ has a node at 0.64 fm. This gives a reasonable overlap with 
the expected wave function for though by no means perfect. This discussion is only semi- 
quantitative because actual radial wave functions are unknown, but provides some explanation 
why the decay to J/^/n may be weak. 

The next task for experiment is a spin analysis of Z(4430), so expressions will be given for 
partial waves. It is highly desirable to use the amplitude for the full process B — > K^'ir; both 
production and decay in principle determine J p , but angular correlations between them provide 
further delicate information which constrains the analysis strongly. My experience is that 100 
fully analysed events achieve the same significance level (in terms of \ 2 °f the fit to data) as 
400-800 events where only information from decays is used. 

Formulae simplify greatly in the rest frame of the Z(4430). There, the angular momentum of 
the kaon in the production reaction is expressed by its 3-momentum K. Likewise, the angular 
momentum of the pion in the decay Z — > tp'n is given by its 3-momentum IT. However, a small 
correction is required for the Lorentz transformation of spin between its rest frame and that 
of Z(4430). Suppose the xz plane is defined by the recoil pion and leptons from the decay of 
Let R be the angle between the (unmeasured) polarisation vector and this plane. In this plane, 
let 9 be the angle between the lepton axis and the pion in the rest frame of . The polarisation 
vector e of the VP' is orthogonal to the lepton axis in the rest frame and may be written as a 
4-vector e = [cos R, sin R, 0, 0] . The intensity is obtained using < cos 2 R >=< sin 2 R >= 1/2. 
A simple trick avoids this integration. The e x and e y components may be replaced by 1/V2 
and i/y/2. Intensities are then obtained by taking the modulus squared of matrix elements. 
Appendix 1 of [35J shows that the Lorentz boost between the rest frames of \&' and Z(4430) 
gives a polarisation vector e' = (1/\^2)[1 + (7 — 1) sin 2 9, i, sin 9 cos#(l — 7)?, (3"f sin 9 cos 9} , where 
(3 and 7 are the usual parameters of the Lorentz transformation. In the Z rest frame, the fourth 
component of e! drops out of all matrix elements. Numerically, (7 — 1) = 0.064, so the effect of 
the Lorentz boost is small. 




(28) 
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After this preliminary, expressions for matrix elements M in the Z rest-frame are simple. 
For J p = _ , M = e'.n and K does not contribute. The choice of xz plane gives U y = 0, so 
the angular distribution is da/dQ oc (e' x Il x ) 2 + (e' z Il z ) 2 ~ n 2 . If the small term in (7 — l) 2 is 
dropped, this is proportional to sin 2 9, as Rosner remarks. 

For J p = 1~, M oc K.e! A II = —e'.K A II, so the kaon and pion are preferentially orthogonal 
in the rest frame of Z. For J p = 2", M oc r^T^, where r a/3 = K a K p - (1/3)(K 2 + K 2 y + K 2 ), 
representing L = 2 for the kaon; T a p = e' a Iip + e^II a — (2/3)(e' x II x + e' z Il z ), remembering that 
Uy = 0. If only one spin is present, ~ 100 events may well be sufficient to identify the spin. 
More than one spin would suggest a non-resonant cusp. 

6 Other cases 

The association of some peaks with thresholds may be numerical accidents. Hence examples dis- 
cussed here are not comprehensive, but concentrate on cases where the agreement is remarkably 
close or is otherwise of special interest. 

A narrow peak is observed in D°p by Babar [3B] and Belle [37] with a mass of 2939.8±1.3±1.0 
MeV and T = 17.5 ± 5.2 ± 5.9 MeV and also in E c (2455)7r. Its mass is just below the D°p 
threshold at 2944.9 ± 0.4 MeV. For an S-wave threshold, J p = (l/2)~ or 3/2". It could be 
a molecular state [38]; alternatively a nearby regular Aq state may be captured by the D°p 
threshold. 

The A£(2595) has a mass of 2595.4±0.6 MeV [33], very close to that of Xj+vr-, namely 2593.6 
MeV. It is too narrow to be a pure cusp, but may be a resonance attracted to the threshold 
where decays to Actt are possible. 

The #(4039) lies close to the D*D* threshold at 4014 MeV. However, it has a width of 80± 10 
MeV, making its association with a cusp effect questionable. Dunynskiy and Voloshin discuss 
in detail the complicated dependence on mass of D*D*, D°D° and D S D S channels [40J. 

The -D s (2315) is very narrow and ~ 50 MeV below the DK threshold. It is therefore unlikely 
to be a molecular state, though it could be an example of a diquark-antidiquark configuration. 
In the absence of other clear examples of strongly bound 4-quark systems, the likely explanation 
is a cs state. 

Amongst the light mesons, the /2(1565) has a phase variation which is well determined by 
Crystal Barrel data on pp —* 3tt° at rest. Those data have been fitted simulaneously with data 
on pp — > UJU7T at rest [41] and definitely require a resonance accurately at the uu threshold. 
The /2(1640) reported by the Particle Data Group is the decay to lulu, which is moved upwards 
from 1565 MeV by ujuj phase space. There is no other candidate for the radial excitation of 
/ 2 (1565), so it appears that this state has been captured by the ujuj threshold. Note, however, 
that the well known / (1500) is not attracted to that threshold; its decays into ujuj are weak. 

The p(1900) of the Particle Data Tables has a rather narrow width. This is suggestive of a 
cusp due to the strong pp threshold. 

The 71"! (1400) could be a bare cusp. It has been assumed by most groups to be a resonance. 
However, it is close to the /i(1285)7r and 6i(1235)7r thresholds, which would appear in S-wave 
decays. Dzierba et al. question whether it is a threshold effect or a resonance [12]. The 7Ti(1600) 
appears at higher mass with the same quantum numbers and is observed dominantly in &i(1235)7r 
and less strongly in /i(1285)7r. So the 7i"i(1400) could be a molecular configuration coupled to 
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these thresholds or could be simply a cusp effect. A full analysis is needed of the two alternatives 
using analytic forms like those presented here. The data on pp — > TT°ir r] show only a very weak 
signal in the t]ir P-wave, insufficient to tell the difference between the two alternatives [43J. 

Valcarce, Vijande and Barnea have made an interesting study of mixing between diquark and 
tetraquark configurations [14] . though they do not specifically take the cusp effect into account. 

6.1 X(1812) 

An intriguing case is the sharp ucj) signal reported by the BES collaboration [35], peaking at 1812 
MeV, just above the <puj threshold at 1801 MeV. If it were purely a threshold effect, it should 
peak considerably higher. It is therefore very likely associated with the /o(1790), a resonance 
clearly separated from /o(1710) in BES data on J/ty — > ujKK [36] and 07T7T [UJ. Data on 
the ujKK channel display a strong / (1710) peak, but nothing is visible in unm, despite large 
statistics. Conversely, the / (1790) appears clearly in (pirn, but any (pKK signal is weak. There 
is a factor 22 difference in decay branching ratios to tttt and KK, hence requiring separate 
/o(1710) and / (1790). The / Q (1790) is also observed in J/V -> 7 47r [37] [3H]. It is accomodated 
naturally as the radial excitation of /o(1370). 

The cj)uo decay can arise naturally from a glueball component in / (1790) [49j [SD]- A glueball 
is a flavour singlet. It has flavour content 

F = (uu + dd + ss) (uu + dd + ss). (29) 

If the decay is to vector mesons, the component (uu + dd) (uu + dd) makes three charge combina- 
tions of pp and one of uouj. The component 2 (uu + dd) ss can make 4u<f) or 2(K*°K*° + K* + K*~) 
or some linear combination. 

There are BES I data on J/ty — > r j(K +r K^ K~n + ) showing that the channel j(K*K*) does 
not contain any significant + signal [51J. The paper says: 'Contributions from ++ and 4 ++ 
are small or absent'. A signal with the same magnitude as that of J/^> — > ^(uj(p) in [35] would 
be rather conspicuous near 1800 MeV, because of the small phase space at that mass for K*K*. 
Its absence there may be qualitatively attributed to the fact that (uu + dd)ss has larger phase 
space in KK decays than K*K*. 

7 Concluding Remarks 

A sharp threshold generates a cusp in the real part of scattering amplitudes at the opening of 
a new threshold. This is a dispersive effect, arising from analyticity. If there is an attractive t- 
or w-channel exchange, it can add coherently to the cusp effect and generate a resonance. The 
/o(980) appears to behave in this way. The cusp can also add coherently to the confinement 
'potential' and attract a regular quark resonance to the threshold. This explains why states like 
/ 2 (1565), ^ (1430) and A c (2940) appear at thresholds. 

At the threshold, zero point energy is minimised by the long-range tail of the wave function. 
Mixing between quark configurations and meson-meson states minimises the energy of the linear 
combination in a way analogous to the formation of a covalent bond in chemistry. 

Diquark-antidiquark resonances may exist. However, if they lie close to a threshold, zero-point 
energy will necessarily mix a large molecular component into the wave function. The /o(980), 
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for example, has a component of ~ 60% KK. Section 2 used known parameters of /o(980) to 
examine how the second-sheet pole position is affected by perturbations to its parameters. The 
conclusion is that the cusp mechanism can attract a resonance over a surprisingly large mass 
interval of order ±100 MeV. If the open channel (tctt) is switched off, the resonance becomes a 
bound-state pole just below threshold. 

The X(3872) is too narrow to be fitted as a pure cusp. Although parameters of the cusp 
can be fine-tuned to fit the line-shape observed in decays to pJ/^, they then fail to reproduce 
the peak observed 3.5 MeV higher in DD* decays. Both peaks may be fitted with a resonance 
or virtual state. A definitive distinction between these possibilities requires phase information 
from interferences in Dalitz plots: for a virtual state, the phase falls above the threshold. The 
natural explanation of X(3872) is that the ss 3 Pi radial excitation has been attracted to the 
DD* threshold. The weak decays of cc to non DD channels like pJ/^, cuJ/^, [xo"]l=o and 
[vcc]l=i lead to a very narrow resonance. 

The Z(4430) can be fitted as a resonance. It is quite possible that meson exchanges generate 
sufficient attraction to turn the cusp into a resonances. However, it cannot presently be excluded 
that a non- resonant cusp fits the data using equations (26) and (27). The availability of many de- 
excitation channels such as [DD*]l =1 necessarily produces a strong cusp. Partial wave formulae 
to assist separation of spin-parity assignments 0~, 1 + and 2~ are given in Section 4. The presence 
of more than one J p would suggest a bare cusp. 

Equations (8), (9) and (14) for a sharp threshold improve on the Flatte formula and are just 
as easy to use. The experimental data for / (980) conform with the line-shape predicted by 
these equations. However, as Tdrnqvist remarks, increasing the precision of the formula and its 
parameters may be an academic exercise. A case which does require study with full inclusion of 
the dispersive effect is 7Ti(1405), which could be a resonance or could be just a threshold cusp. 
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